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Abstract 


Application of Asymptotic Waveform Evaluation (AWE) is presented in conjunction with 
a hybrid Finite Element Method (FEM) / Method of Moments (MoM) technique to calculate the 
input characteristics of cavity-backed aperture antennas over a frequency range. The hybrid FEM/ 
MoM technique is used to form an integro-partial -differential equation to compute the electric 
field distribution of the cavity-backed aperture antenna. The electric field, thus obtained, is 
expanded in a Taylor series around the frequency of interest. The coefficients of Taylor series 
(called “moments”) are obtained using the frequency derivatives of the integro-partial-differential 
equation formed by the hybrid FEM/MoM technique. Using the moments the electric field in the 
cavity is obtained over a frequency range. Using the electric field at different frequencies, the 
input characteristics of the antenna are obtained over a wide frequency band. Numerical results 
for an open coaxial line, probe fed cavity, and cavity-backed microstrip patch antennas are pre- 
sented. Good agreement between AWE and the exact solution over the frequency range is 
observed. 
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List of Symbols 


v 

V' 


Del operator 

Del operator over the source coordinates 
Dielectric permittivity of the medium in the cavity 
Dielectric permittivity of the medium in the coaxial feed line 
Kronecker delta defined in equation (22) 

Dielectric permeability of the medium in the cavity 
p -coordinate of the cylindrical coordinate system 
Unit normal vector along the p -axis 
Angular frequency 

Asymptotic Waveform Evaluation 

d q 

qih derivative of A (k) with respect to k ; A ( k ) , evaluated at k 

dk q 

Reflection coefficient at the input plane S inp 

Excitation vector 

d q 

q\h derivative of b (k) with respect to k ; b ( k ) , evaluated at k o 

dk q 

Surface integration with respect to observation coordinates 
Surface integration with respect to source coordinates 
Electric field 


E inp Electric field at the input plane S { 

e (k) Electric field coefficient vector 

e . Incident electric field due to coaxial line at the surface S 

inc inp 

e , Reflected electric field into the coaxial line at the surface S 

ref inp 

H fl/j Magnetic field at the surface S ap 
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Magnetic field at the surface S jnp 
Frequency 

Wavenumber at any frequency / 

Wavenumber at frequency f 
Magnetic current at the surface S ap 

«th moment of AWE ( n=0,],2,3,4 ) 

Normal unit vector 
Factorial of number q 

Distance between the source point and the observation point 

Radius of inner conductor of the coaxial feed line 

Radius of outer conductor of the coaxial feed line 

Vector testing function 

Vector testing function at the surface S ap 

Normalized input admittance of the antenna 


Unit normal along Z-axis 



1. Introduction 


Cavity-backed aperture antennas are very popular in aerospace applications due to their 
conformal nature. These antennas can be analyzed using the integral equation or differential 
equation methods. The integral equation approach involves the solution of a fully dense matrix 
equation and mathematically complex for inhomogeneous material and arbitrarily shaped 
cavities. The differential equation method can easily handle the arbitrarily shaped cavities with 
inhomogenous materials, but requires boundary truncation. Hybrid techniques have become 
attractive for numerical analysis of these type of problems due to their ability to handle arbitrary 
shape of the cavity and complex materials that may be required for the antenna design. The 
combined Finite Element Method (FEM) and Method of Moments (MoM) technique in particular 
has been used to analyze various cavity-backed aperture antennas[l,2]. In the combined 
FEM/MoM technique, FEM is used in the cavity volume to compute the electric field, whereas 
MoM is used to compute the magnetic current at the aperture. Using Galerkin’s technique and 
forming simultaneous equations, the electric field is solved. For the combined FEM/MoM 
technique, the cavity is divided into tetrahedral elements and the aperture is discretized by 
triangles. Simultaneous equations are generated over the subdomains and are added to form a 
global matrix equation. This results in a partly sparse and partly dense symmetric complex matrix, 
which can be solved either by a direct solver or by an iterative solver. The electric field hence 
obtained is used to compute the radiation characteristics and input characteristics of the antenna. 

In most practical applications, input characteristics such as input impedance or input 
admittance are of interest over a frequency range. To obtain the frequency response of the 
antenna, one has to repeat the above calculations at every incremental frequency over the 
frequency band of interest. If the antenna is highly frequency dependent, one needs to do the 
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calculations at fine increments of frequency to get an accurate representation of the frequency 
response. This can be computationally intensive for electrically large cavity with a large aperture 
and in some cases computationally prohibitive. To alleviate the above problems, the application 
of Asymptotic Waveform Evaluation (AWE) [3] to the combined FEM/MoM technique is 
proposed. Recently a detailed description of AWE was applied to frequency domain 
electromagnetic analysis and is presented in [4], AWE has also been used to predict radar cross 
section (RCS) of perfect electric conductor (PEC) bodies over a frequency range [5], 

In this report, we describe the application of AWE for calculating the input characteristics 
of a cavity-backed aperture antenna over a band of frequencies using the combined FEM/MoM 
technique. In the AWE technique, the electric field is expanded in Taylor’s series around a 
frequency. The coefficients of Taylor series (called ’moments’) are evaluated using the frequency 
derivatives of the combined FEM/MoM equation. Once the moments are obtained, the electric 
field distribution in the cavity can be obtained at any frequency over the frequency range. Using 
this field distribution, the input characteristics of the cavity-backed aperture antenna are obtained. 

The rest of the report is organized as described below. In section 2, the combined 
FEM/MoM formulation is presented. In section 3, AWE implementation for the combined 
FEM/MoM technique is described. Numerical results for an open coaxial line, a coaxial cavity, 
and a cavity-backed microstrip patch antenna are presented in section 4. The numerical data are 
compared with the exact solution over the bandwidth. CPU time and storage requirements for 
AWE formulation are given for each example and are compared with those required for exact 
solution at each frequency. Concluding remarks on the advantages and disadvantages of the AWE 
technique are presented in section 5. 
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2. Combined FEM/MoM Technique for Cavity-Backed 
Aperture Antennas in Infinite Ground Plane 

The geometry of the problem to be analyzed is shown in figure 1. For linear, isotropic, and 
source free region; the electric field satisfies the vector wave equation: 

Vxf— VxE^ -& 2 £ r E = 0 (1) 

where p, r , e r are the relative permeability and relative permittivity of the medium in the cavity. 
The time variation exp (/to/) is assumed and suppressed throughout this report. Applying the 
Galerkin’s technique, equation (1) can be written in “weak form” as [6] 

JJj (VxT) .('ivxEW* 2 E r JJjT.Et/v-yconJ j (Txn) -H ap ds 

v v^ r y v 5 

“p 


= /conj JT. (nxH inp ) ds (2) 

np 

where T is the vector testing function, S ap is the aperture surface, and S jnp is the input surface 
(see figure 1). is the magnetic field at the aperture and H inp is the magnetic field at the input 
surface. 

In accordance with the equivalence principle [7], the fields inside the cavity can be 
decoupled to the fields outside the cavity by closing the aperture with a PEC and introducing the 
equivalent magnetic current. 

M = E x z (3) 

over the extent of the aperture. Making use of the image theory, the integrals over S ap in equation 
(2) can be written as 
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M*J J (T x n) • H ap ds 


s/Jv 


exp <-m ds . 




V 


ds 


J 


rJj <v * T -> )JJ (V'.M)2Efca*u.j* 


( 4 ) 


where T s = Txn and R is the distance between source point and the observation point. V' indi- 
cates del operation over the source coordinates and ds' indicates the surface integration over the 
source region. 

Though the analysis presented in this report is not restricted to any specific input feed 
structure, we restrict the presentation of the formulation to the coaxial line as the input feed 
structure. The cross section of the coaxial line is shown in figure 2. Assuming that the incident 
electric field is the dominant transverse electric and magnetic (TEM) mode and the reflected field 
also consists of TEM mode only, the electric field at the input plane 5 f . is given by 

E inp = e inc eX P[-J k fi^c Z ) + e ref eX p[j k J^c Z ) ( 5 ) 

where 


and 


e = p- 

inc r 


271 In 


v r iy 


( 6 ) 


'ref 


= a e 
o inc 


a o is the reflection coefficient and is given by 


( 7 ) 
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r 2 is the outer radius and r, is the inner radius of the coaxial line. e rc is the relative permittivity 
of the coaxial line. 

Using equation (5) to calculate H inp , the surface integral over S inp in equation (2) can be 
written as 

Wj | T * (" XH inp) ds 

S,„„ 



Substituting equation (4) and (8) in equation (2), the system equations for the combined FEM/ 
MoM technique can be written as 

III if (VxT) • ( VxE) dv-/c 2 £ r JJ|T • E<iv 

v Vr y 
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The volume of the cavity is subdivided into small volume tetrahedral elements. The 
electric field is expressed in terms of the edge vector basis functions [6], which enforce the 
divergenceless condition of the electric field explicitly. The vector testing function is also 
expressed in terms of the edge vector basis functions following the Galerkin’s method. The 
discretization of the cavity volume into tetrahedral elements automatically results in discretization 
of the surfaces S ap and S jnp into triangular elements. The volume and surface integrals in 
equation (10) are carried out over each element to form element matrices and the element 
matrices are assembled to form global matrices. Equation (10) can be written in matrix form as 

A (k) e (k) = b(k) (11) 

A ( k ) is a partly sparse, partly dense complex symmetric matrix, b(k) is the excitation vector, and 
e(k) is the unknown electric field coefficient vector. A(k) is evaluated as a sum of three matrices. 

A(k) = A ; (k) +A 2 (k) +A 3 (k) +A 4 (k) (12) 

where 


A x (k) = jjj — (VxT) • ( VxE) dv-k 2 e r ^jjj • Edv 


■V*> = 


S ap V 


^3 (*) = yJ J ( ] J J (V'.M) exp ds ' 


(13) 

(14) 

(15) 
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The matrix equation (11) is solved at any specific frequency, f g (with wavenumber k g ) 
either by a direct method or by an iterative method. The solution of the equation (11) gives the 
unknown electric field coefficients which are used to obtain the electric field distribution. Once 
the electric field distribution is known, the input reflection coefficient can be calculated using 
equation (8). The input plane is placed at = 0 , and the reflection coefficient is calculated as 



The normalized input admittance at S- is given by 


Y = LzL 
“ i + r 


The input admittance given in equation (18) is calculated at one frequency. If one needs the input 
admittance over a frequency range, this calculation is to be repeated at different frequency values. 


3. AWE Implementation 

The general implementation of AWE for any frequency domain technique used for 
electromagnetic analysis is given in detail in [4]. As shown in the previous section, the solution of 
equation (11) gives the unknown electric field coefficient vector e ( k o ) at a particular frequency 
f o . However at any k, e (k) can be expanded in Taylor series as 
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oo 


( 20 ) 


e <*> = 


n = 0 


with the moments M n given by [4] 


-l 


M=A <*„) 


b"" (k ) " (1-8 )A (,> (k)M 

x o' _ ^ go’ v o ' n - 1 

n\ ~ ^ 

<7 = 0 


q- 


( 21 ) 


(<?) 


A ( k o ) is the ^th derivative with respect to k of A(k) given in equation (12) and evaluated at 


(9) 


k Q ■ Similarly, b (k o ) is the qth derivative with respect to k of b(k) given in equation (16) and 


evaluated at £ . The Kronecker delta 6 is defined as 

o go 


8 = { 

go L 


1 q = 0 

0 q £ 0 


( 22 ) 


The q th derivatives of A(k) and b(k) are evaluated and are given in detail in the Appendix. 

Once the moments of AWE are obtained, the electric field coefficients at frequencies 
around the expansion frequency are obtained by using the equation (20). The electric field hence 
obtained is used to compute the input characteristics of the cavity-backed aperture antenna over a 
frequency range. 


4. Numerical Results 

To validate the analysis presented in the previous sections, a few numerical examples are 
considered. Input characteristic calculations over a frequency range are done for an open coaxial 
line, coaxial cavity, and cavity-backed square and circular microstrip patch antennas. The 
numerical data obtained using AWE are compared with the results calculated at each frequency 
using the computer code CBS3DR[9], which implements the combined FEM/MoM techniquef2]. 
We will refer to the latter method as “exact solution”. From section 3, it can be observed that the 
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inverse of matrix A ( k o ) is found once and is used repeatedly to find AWE moments. Due to the 
hybrid FEM/MoM technique, matrix A ( k o ) is partly sparse and partly dense. The Complex 
Vector Sparse Solver (CVSS) [10] is used to LU factor the matrix A ( k Q ) once and the moments 
are obtained by backsolving the equation (21) with multiple righthand sides. All the computations 
reported below are done on a CONVEX C-220 computer. 

(a) Open Coaxial line: 

An open coaixial line radiating into an infinite ground plane (fig. 3a) is considered. A 
finite length of the line is used for FEM discretization. The input plane S ( is placed at z = 0 
plane and the radiating aperture at z = 1 cm . The discretization of the coaxial line resulted in 
1119 total unknowns, and the order of the dense matrix due to MoM is 144. The frequency 
response of the input admittance is calculated with 6 GHz as the expansion frequency. The AWE 
moments are calculated at 6 GHz and are used in the Taylor series expansion. The frequency 
response from 4 GHz to 8 GHz is plotted in figure 3(b) along with the exact solution calculated at 
different frequencies. A good trend in frequency response is predicted by the 2nd order AWE 1 , 
whereas a very good agreement can be seen between the 5th order AWE frequency response and 
the exact solution over the frequency range. The 1119X1119 hybrid FEM/MoM matrix exact 
solution took around 150.7 secs of CPU time to fill the matrix and 10 secs to LU factor the matrix 
at each frequency. The 5th order AWE frequency response calculation took 328 secs of CPU time 
to fill the matrices including the frequency derivative matrices and 10 secs to LU factor the 
A (k o ) matrix. The exact solution was carried out at nine frequency points with (160.7X9) 
1446.30 secs of total CPU time. With AWE, the frequency response was calculated with 0.1 GHz 
frequency increments. It can be seen that there is a substantial amount of savings in CPU time by 

1. As AWE is a purely mathematical approximation to the solution, it is observed that at some frequencies, 

2nd order AWE results in unrealistic values of conductance. 
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using AWE, when frequency response of input characteristics are required with fine frequency 
increments. 

(b) Open Coaxial Cavity: 

An open coaxial cavity fed by a 50Q coaxial line (fig. 4) is considered as a second 
example. The input plane S inp is placed at z = 0 plane and the radiating aperture at 
Z = 0.952 cm plane. The cavity volume is discretized using tetrahedral elements, which resulted 
in 4541 total unknowns and the order of the dense matrix due to MoM is 666. The frequency 
response of the return loss (=201og|T| ) is calculated with 6 GHz as the expansion frequency. The 
AWE moments are calculated at 6 GHz and are used in the Taylor series expansion. The 
frequency response from 5 GHz to 7 GHz is plotted in figure 5 along with the exact solution 
calculated at different frequencies. It can be seen from figure 5 that 2nd order AWE could not 
predict the frequency response over the frequency range, whereas a very good agreement can be 
seen between the 5th order AWE frequency response and the exact solution over the frequency 
range. The 4541X4541 hybrid FEM/MoM matrix exact solution took around 2027 secs of CPU 
time to fill the matrix and 264 secs to LU factor the matrix at each frequency. The 5th order AWE 
frequency response calculation took 4867 secs of CPU time to fill the matrices including the 
frequency derivative matrices and 264 secs to LU factor the A (kj matrix. The exact solution 
was carried out at nine frequency points with (2,291X9) 20,619 secs of total CPU time. With 
AWE the frequency response was calculated with 0.1 GHz frequency increments with 5,140 secs 
of total CPU time. 

(c) Cavity-Backed Square Microstrip Patch Antenna: 

A cavity-backed square microstrip antenna radiating into an infinite ground plane (fig. 6) 

is considered. The input plane S. is placed at z = 0 plane and the radiating aperture at 
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Z - 0.16cm. The discretization of the cavity volume resulted in 2,160 total unknowns and the 
order of the dense matrix due to MoM is 544. The frequency response of the input impedance 
( 1 / T. ) is calculated with 4 GHz as the expansion frequency. The AWE moments are calculated 
at 4GHz and are used in the Taylor series expansion. To obtain an accurate frequency response 
over a wider frequency range, another set of moments are calculated at 4.3 GHz. The frequency 
response from 3.8 GHz to 4.5 GHz is plotted in figure 7 along with the exact solution calculated at 
different frequencies. A very good agreement can be seen between the 5th order AWE frequency 
response and the exact solution over the frequency range. The 2160X2160 hybrid FEM/MoM 
matrix exact solution took around 1558secs of CPU time to fill the matrix and 122secs to LU 
factor the matrix at each frequency. The 5th order AWE frequency response calculation took 3754 
secs of CPU time to fill the matrices including the frequency derivative matrices and 122 secs to 
LU factor the A ( k o ) matrix. The exact solution was carried out at nine frequency points with 
(1,680X9) 15,120 secs of total CPU time. With two expansion points, AWE took 7,752s ecs of 
total CPU time. With AWE the frequency response was calculated with 0.01 GHz frequency 
increments. 

(d) Cavity-Backed Circular Microstrip Patch Antenna: 

A cavity-backed circular microstrip antenna radiating into an infinite ground plane is 
shown in figure 8. The input plane S- is placed at z = 0 plane and the radiating aperture at 
z = 0.16cm. The discretization of the cavity volume resulted in 6,363 total unknowns and the 
order of the dense matrix due to MoM is 469. The frequency response of the input impedance 
( 1/y. ) is calculated with 6 GHz as the expansion frequency. The AWE moments are calculated 
at 6 GHz and are used in the Taylor series expansion. To obtain an accurate frequency response 
over a wider frequency range, another set of moments are calculated at 5.8 GHz. The frequency 
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response from 5.6 GHz to 6.2 GHz is plotted in figure 9 along with the exact solution calculated at 
different frequencies. A very good agreement can be seen between the 5th order AWE frequency 
response and the exact solution over the frequency range. The 6,363X6,363 hybrid FEM/MoM 
matrix exact solution took around 1,250 secs of CPU time to fill the matrix and 112secs to LU 
factor the matrix at each frequency. The 5th order AWE frequency response calculation took 
2,967 secs of CPU time to fill the matrices including the frequency derivative matrices and 112 
secs to LU factor the A ( k o ) matrix. The exact solution was carried out at seven frequency points 
with (1362X7) 9,534 secs of total CPU time. With two expansion points, AWE took 6,178secs of 
total CPU time. With AWE the frequency response was calculated with 0.01 GHz frequency 
increments. Considering the fact that with AWE around 60 frequency calculations could be 
carried out with less CPU time compared to calculate 7 frequency points with exact solution, 
AWE has a distinct advantage and is essential if one has to determine the exact resonant 
frequency. 


5. Concluding Remarks 

The AWE technique is applied to the hybrid FEM/MoM technique to obtain the frequency 
response of the input characteristics of cavity-backed aperture antennas. The frequency response 
of input characteristics of an open coaxial line, coaxial cavity, square microstrip patch anetnna, 
and a circular patch antenna are computed and compared with the exact solution. From the 
numerical examples presented in this work, AWE is found to be superior in terms of CPU time to 
obtain a frequency response. It may be noted that although calculations are done in frequency 
increments of 0.1 GHz or 0.01 GHz for the examples presented, the frequency response at even 
finer freqeuncy increments can also be calculated with a very nominal cost. The application of 
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AWE to three dimensional cavity-backed aperture antennas (without the infinite ground plane) is 
of interest for future research. The accuracy of AWE over a desired frequency band and its 
relation to the order of AWE to be used are also of interest for future research. With all of these 
topics addressed, AWE will be a good computing tool for the design of cavity-backed aperture 
antennas. 
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Appendix 

Derivatives of A(k) and b(k) w.r.t. k 


The frequency derivatives of A(k) and b(k) are evaluated and are given below. From 
equation (12): 


A (q) (k) = 

dk q 

= A (q) (k) +A ( 2 q) ( k ) +A\ q) (k) + A {q) (*) q=0,l,2,3, 

(A.l) 

From equation (13) 



A, (0) (k) 

= JJJ — (VxT) • (VxE) Jv-£ 2 e r JJjT • Edv 

V y 

(A. 2) 


A, (,) (£) = -2*eJJJt«E</v 

V 

(A. 3) 


a\ 2} (k) = -2e r JjjT • Edv 

V 

(A.4) 


S 

w 

II 

o 

*0 

IV 

U> 

(A.5) 


From equation (14) 


V 

+ (A.7) 

S ap \ S ap J 

j 1 M ( 2n) [ 2) ! ( ~ jR) q 3 + 2 ?* (-•/*)’ 2 + k 2 (-jR) q 'j exp (-jkR) ds')ds for q> 1 ( A. 8) 
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From equation (16) 


A 3 0) (k) = h.1 J ( V#T ^ ) 1 1 J ( V ' • M) exp ( R jkR) ds' 1 ds (A.9) 


A\ q) (k) = JJ (V.T 5 ) { J J (V' • M) ( J (~jR) (q ~ 1} ex P (~jkR) ds'} ds (A.10) 


From equation (17) 


■-7§-{// T *(pW{iJ E *(p) b 

2nln — p ^ 5 .» P J ^ s in P 


(A. 11) 


a: ; (k) = 


a; ; (k) 


(A. 12) 


Ar (k) = 0 


q>2 


(A. 13) 


From equation (18) 


]i rc pita — s “' 

V l r i J 


(A. 14) 


* W <*) = (- 2^1 JT 1 -^—’ l fe(0) (*) 

V ' L Jtftrc* J 


(A. 15) 


Equation (A. 15) is written in a compact form, however, it must be simplfied before evaluating at 


z, = 0. 
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Figure 1 Geometry of a cavity backed aperture in finite ground plane. 
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Input 

(a) 



Figure 3 (a) Open coaxial line in an infinite ground plane. Inner radius Outer 

radius r 2 = 1.57cm, e r =1.0 and L=1.0cm 
(b) Normalized input admittance as a function of frequency. 
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Figure 4 Geometry of a coaxial cavity in an infinite ground plane. Outer radius of the 
coaxial cavity=l”, Inner radius of the coaxial cavity=0.0181” and L=3/8”. The 
cavity is fed by a 50£2 coaxial line. 
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Frequency (GHz) 


Return loss versus frequency of the coaxial cavity (figure 4). 





Y 

A 



i_i 

50Q coaxial feed 


Figure 6 Cavity-backed square microstrip patch antenna in an infinite ground plane fed by a 
50£2 coaxial line. 
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Normalized input impedance versus frequency of the cavity-backed square 
microstrip antenna (figure 6). 







Figure 8 



Cavity-backed circular microstrip patch antenna in an infinite ground plane fed by 
a 50Q coaxial line. 
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Figure 9 Normalized input impedance versus frequency of the cavity-backed circular 
microstrip antenna (figure 8). 
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